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ABSTRACT 


Turbulence in the atmosphere drives the formation of temperature 
inhomogeneities that scatter and diffract propagating electromagnetic waves, adversely 
affecting laser weapons and high-resolution optical systems. Military operations require 
reliable turbulence profiles for the development and validation of turbulence prediction 
models.. 

This research investigated the false turbulence contribution caused by well-known 
temperature steps in the vertical profile of the atmosphere, especially in the stratosphere. 
The homogeneity and isotropy requirements of structure functions were used to develop a 
technique to remove the false contribution to the temperature structure constant, C,”. 

Both 1.54 cm and 5.82 m saad: ssschition profiles with 0.001 to 0.01 K 
temperature resolution were collected from a balloon flight. Steps of 0.1 to 1 K in the - 
vertical temperature profile produce abrupt changes in the mean temperature that obscure 
the measurement of the actual turbulent fluctuations. Removing these anomalies exposed 
the underlying C,’ distribution. Application of the new technique for several sampling 
intervals revealed a Kolmogorov inertial subrange extending from ~25 cm to ~10 m. The 


potential of this technique to compute the isoplanatic angle, 6, coherence length, 7,, and 
Greenwood frequency, f,, reliably by using inexpensive balloons should benefit airborne 


and space-based laser programs. 
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I. INTRODUCTION 


Atmospheric turbulence degrades propagating electromagnetic waves by the 
processes of scattering and diffraction. Many laser and telescope systems of interest to 
the military operate within the visible through near infrared portion of the spectrum, 


referred to as the optical spectrum. Systems that require less than 1-10 firad angular 


resolution are especially prone to degradation from atmospheric turbulence over paths of 
~1 km or greater. | 

The Airborne Laser (ABL), sponsored by the USAF, is a system that transmits a 
high-energy beam over several hundred kilometers and focuses it to a small — size. 
The effects of atmospheric turbulence are a primary limitation to the system performance. 
Measuring the distortion of the wavefront in near real time and compensating with a 
deformable mirror can increase the energy density to the level needed to destroy a missile 
during the boost phase. 

The ability to forecast the severity of atmospheric turbulence accurately would be 
of value in mission planning and coordinating real time operations for ABL and other 
military adaptive optics systems. Developing a reliable forecast model requires accurate 
validation data for comparison. The potential of inexpensive meteorological balloons to 
provide accurate turbulence profiles from the temperature profile is enticing. 

A VIZ microsonde package launched in March 1995 from Wichita, Kansas, was 


augmented with two additional microthermal temperature instruments. They provided 


1.54 cm and 5.82 m vertical resolution profiles of temperature in addition to the standard 








pressure, temperature, humidity, wind speed and wind direction measurements. The 
weather consisted of dry snow at low altitudes topped by clear air. 

There was cross coupling between the VIZ microsonde electronics package and 
the high-resolution thermocouple electronics at the 1.25 sec sample period that produced 
anomalous temperature spikes. These were removed by data processing and editing. 
Blind application of a second-order structure function to the temperature profiles revealed 
the signature of temperature steps and confirmed the tendency of the stable stratosphere 
to organize itself into a stepped profile of temperature on all scales. The homogeneity 
and isotropy requirements of structure functions led to the development of a quantitative 
measure and procedure to remove the false turbulence caused by the temperature steps. 
Confirmation of this technique came from the evidence of a Kolmogorov type inertial 
subrange after removal of the temperature steps. The inertial subrange extended from 
~25 cm to ~10 m. This research produced a reference turbulence profile and 
measurement procedure suitable for validation of forecast models. Horizontally spaced 
probes are also subject to these same step effects. Further research may show that a 
correction factor can remove these effects, or that smaller horizontal separations may be 
necessary. The opportunity exists for ABL and other programs to take advantage of the 
narrow Kolmogorov window with inexpensive 5 m vertical resolution instruments. 

Chapter II presents a background of theory on optical turbulence, structure 
functions, and optical parameters. The new theory of temperature steps is also presented. 


- The balloon instrumentation and preliminary data processing are discussed in Chapter III. 





Results of the investigation are presented in Chapter IV. Finally, concluding remarks are 


made in Chapter V. 








THIS PAGE INTENTIONALLY LEFT BLANK 








Ii. BACKGROUND 


A. OPTICAL TURBULENCE 


Variations in the index of refraction of a medium distort the wavefront of 
propagating electromagnetic waves. Consequently, the constant phase surface of the 
wave departs from a flat surface. Figure 1 shows an incident plane wave passing through 
a region of fluctuating, optically turbulent, refractive index. The wavefront speeds 
up/slows down through lower/higher index of refraction regions. 

Treating the atmosphere as an ideal gas, the index of refraction, n, becomes a 


function of pressure P, temperature T, and humidity Q, 


n=n(P,T,OQ). (1) 
For air, there is also a weak dependence on wavelength, which is negligible over the 
optical spectrum (Jones, 1981). For optical wavelengths, the index of refraction may be 


written as (Tatarski, 1961, pp. 55) 


79x10° 
+ —___—-—— 
T 


n=l ae Z)4 


Fluctuations in P (mbar) and T (K) from the mean values will result in corresponding 
fluctuations in the index of refraction. The incremental change in the index of refraction 


can be found by taking the total differential of n, 


dn -19x10"( apa, (3) 
ar = 





= 


Figure 1. Variations in the index of refraction distort a plane wave as it passes through an 
optically turbulent medium. 








where the two increments represent turbulent departures from the mean. The humidity 
contribution to Equation (2), not shown, has a negligible contribution to Equation (3) and 
was ignored for simplicity. A further simplification arises because the pressure variation 
term is very small compared to the temperature variation term. This is because typical 
atmospheric wind velocities are small compared to the speed of sound, allowing time for 
equilibration, which suppresses the formation of acoustic waves. 

Positive and negative changes, dn, from the mean index of refraction distort the 
phase of a propagating wave, By squaring a simplified form of Equation (3), we form a 
statistical quantity, 


2 
dn® = [79 x10°° z) a (4) 


2 


which is closely related to the variance of the fluctuations. The theory of random 


functions is appropriate in studying the spatially or temporally averaged temperature 
term, (aT *) , and is discussed in the next section. 
B. | STRUCTURE FUNCTIONS 


1. First Order 


The variance of a stationary variable, such as temperature, can be written as 


(dT?) = ((r-(r))). (5) 
One difficulty with this approach is that the average temperature may vary over the 


region of interest or over time; it is not stationary. When this occurs, the variance 











includes a contribution from the overall change in mean temperature and no longer 
represents only temperature fluctuations. 

Structure functions provide a more robust method to examine the temperature 
fluctuations than the variance. The 1° order temperature structure function is a function 


of the temperature at two positions, 1 and 2, and is written as 


Dr =((7,-1,)), (6) 
where the superscript indicates the function is a first order difference. It quantifies the 
level of temperature fluctuations without requiring stationarity. The difference 
completely rejects constant values in the variable, as does the variance in Equation (5). 
The advantage of the structure function in Equation (6) over the variance is that it rejects 
gradual shifts in the mean value of temperature common to both terms. As the separation 
between 7, and 7; increases, the difference in the mean values of 7, and 7, approaches 
and eventually dominates the random fluctuations in temperature. At this distance and 
beyond, the structure function no longer rejects mean value gradients. The steeper the 
gradient in the mean value, the shorter the separation distance over which the structure 
function is useful. Higher order structure functions provide a solution to the mean 
gradient problem and are discussed in Section 2. 

Within regions of limited size, when the turbulence is homogeneous and isotropic 
the structure function of the atmosphere has a Kolmogorov dependence on distance r 


between the points, 1 and 2, given by (Tatarski, 1961, pp. 46) 


DY; (r)=C,77r?”. (7) 


Since this function increases with Separation distance, it is interpreted as a measure of the 
temperature turbulence. For homogeneous turbulence, the structure function is related to 


the correlation function in the following manner, 


D";(r)= 2B, (0)- B, (r)] (8) 
where B, (r) is the correlation function of temperature (Tatarski, 1961, pp. 10). Equation 
(8) provides additional insight into why structure functions begin at zero and rise with r 
monotonically. | 

The index of refraction structure function has the same form as the temperature 


structure function in Equation (7), 


Dp”, (r) os Co! (9) 
where the index of refraction structure constant is related to the temperature structure 


constant by, refer to Equation (4), (Tatarski, 1961, pp. 79) 


n 


C,’ =C,'| 79 tos 2) 10 
=C, Xx T? . (10) 


Zz Higher Order 


Generalized structure functions of n™ order have the advantage that they remove 
polynomial trends in the data of order n-1. This property makes them useful in analyzing 
geophysical variables such as temperature and wind speed. For example, a second-order 


structure function removes both a constant and a linear trend in the data. The second- 


order structure function uses the difference in the slopes of the temperature profile at two 














points, A and B. The slopes at A and B are found from a finite difference in the 


temperature, 


aT | aT 
Ban Ben 


The second-order structure function becomes 


Dd, =([(7, -7,)-(,-T)P) = -24 +00) - (12) 
In general, the n” order structure function (Yaglom, 1987, pp. 427) can be written in 


terms of the function at n+1 different points and the binomial coefficient, 


Dy = (5 (- »(J5]) | (13) 


The theory has an interesting and useful result for structure functions of the form 


Dr(r)=CMr?r", (14) 
under which the Kolmogorov expression in Equation (7) falls. For m>0, the various 


order coefficients obey the relation, 


n-1 ( ‘) 
cl),? cof a "(ns)", (15) 
j=0 


to the first order coefficient (Walters, 1995). The 1" order coefficient is the desired 


quantity, with C,° =C ()? and is found in Chapter IV from the 2™ order structure 


10 





function by the ratio C,” =C");?/2.412. Section D will show the utility of the 


temperature and index of refraction constants. 


3. Inertial Subrange and Power Spectrum 
This section discusses the limits and behavior of the Kolmogorov result in 


Equation (7). The Kolmogorov structure function behavior holds for homogeneous and 
isotropic volumes. Beyond the outer scale length, L,, the assumptions do not hold and 
the law fails. For example, the turbulence in a thundercloud is quite different from the 
quiescent air just outside the boundary of the cloud, introducing anisotropy. As the 
turbulent cascade proceeds to smaller eddies, the viscous dissipation of the air increases 
and ideally prevents turbulent motion below an inner length scale, /,, having converted 
the turbulent kinetic energy of motion into thermal energy of the air molecules. Using a 
power series expansion, the structure function has an r*> power law below /, (Tatarski, 
1961, pp. 32). Figure 2 illustrates the idealized behavior of experimentally determined 
temperature structure functions against the theoretical functions. 

The power spectrum of a turbulent variable reveals the range over which the 
Kolmogorov power law applies. Figure 3 is a representative log-log plot of an idealized 
one-dimensional temperature power spectrum. The inertial subrange is the region with 
constant slope, where the —5/3 slope follows from the r?? form of the structure function 
law (Tatarski, 1961, pp. 25). This result hails from the general equation relating the 


structure function of Equation (14) to the power spectral density W(o), 








-" 
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Figure 2. 


Idealized structure function behavior with sampling interval r. The 
Kolmogorov 2/3 law holds between /, and L,. Quadratic behavior is observed below 
/,. After Tatarski, 1961, pp. 33. 
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Figure 3. Idealized power spectrum of temperature as a function of frequency. The 


inertial subrange, viscous-convective, and dissipation regions are shown (Tennekes, 19, 
pp. 284). 
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W(@)= TR sn Cal, (16) 


for 0<m<2 (Tatarski, 1961, pp. 13). Inner and outer scale lengths are obtained from 
the spatial frequencies that bound the power law region of the spectrum. The dissipation, 
or viscous convective, kink at the right end of the spectrum is a result of the mechanical 
viscosity stopping the turbulent motion before temperature diffusion can equalize the 
temperature. Beyond the kink, the spectrum falls off rapidly as thermal diffusion 


smoothes the variations. 


C. TEMPERATURE STEPS 


The structure function approach can fail when applied over a region that is not 
homogeneous and isotropic. A false contribution to the second-order structure function 
can occur when the gradient of the temperature profile changes appreciably over the 
measurement distance with respect to the fluctuations. Figure 4 presents a simple model 


in which the temperature increases with slope b, below, and slope b, above, height z,. 


The abrupt change in slope alters the behavior of the structure function. 


Consider the temperature, T(z,t), as a function of height, z, and time, r, 


T(z,t)=T(z)+é(2,2), (17) 
to be a sum of the mean, time averaged, temperature, em (z,t)) = T(z), and a fluctuation 


term, E(z, t), representing the turbulence. By this definition, the time average of the 


fluctuation term is zero, 
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To T, . T, T 


Figure 4. A simple model of a temperature step. Second-order structure functions 
computed at the step, slope discontinuity, are artificially high. 


is: 





(E(z,1)), =0. (18) 
Following Figure 4, the linear function of temperature with height is 
= a+b2z, z2=z 
Fizi= : ma ° 19 
(2) ee _ oa 
where a,,a,,5,, and 5, are constants. Continuity of temperature across the step results 


in the boundary condition, 


a, +6,z,=a,+b,z,. (20) 
Substituting the expressions for temperature in Equation (19) into the second-order 


structure function of Equation (12) yields 


Dr(r)= (7, -27, +P) = (I, -b y+ - 28, +&)P)- (21) 
The difference operation eliminates the constants, a, and a,, and the sampling interval 
r=Z,—2Z, =2Z,—Z, = Az is the height increment. Equation (21) now includes terms for 
the slopes and the temperature fluctuations. Rewriting B=b, —b, and multiplying out 


the terms results in, 


Dr (r)= ([BrP }+([(E -28, +4)P }, (22) 
after using Equation (18) to eliminate the cross terms. The second term of Equation (22) 


is the desired 2" order structure function. Using the Kolmogorov expression in Equation 


(7) and the ratio C,” = C7? /2.412, from Equation (15), gives 
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Dr (r) = B’r? +2.41C,7723, (23) 
The measured 2™ order structure function is a combination of the true temperature 
turbulence, obeying a 2/3 law, and a false quadratic slope discontinuity term. A measure 


of the relative contribution by the slope discontinuity term is the ratio 


2 
BY ass 


=r 
2.41C,,” 


(24) 


The second-order structure function must only be calculated within constant gradient or 
homogenous regions under isotropic conditions in order to preclude a false contribution. 


Removing the steps from the structure function should reveal the true turbulence 
constant, C,’. Although the development pertains to 2™ order structure functions, the 
theory is general to n™ order. Slow trends in the mean value may be ignored = iP. 
order structure functions if the corresponding Q factor is much smaller than 1, as was 


mentioned in Section B. Similarly, rapid changes in curvature would add a false 


contribution to 3" order structure functions, etc. 


D. OPTICAL PARAMETERS 
1. Coherence Length 


The coherence length, 7,, or Fried parameter, is useful because it is the equivalent 


maximum aperture an optical system could have and still obtain a near perfect diffraction 


limited image, Pree =Alry vice Q,,,=A/D. Equivalently, the maximum sub-aperture 


size of an adaptive optics system, d , is dictated by the lowest value of 7, for which it is 
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necessary to compensate. Physically, it is a characteristic scale over which the 


electromagnetic field has approximately constant phase. For a plane wave, 


| v7 -3/5 
r, = 2,{221 k? fc." (25) 
0 


where the integral is of the phase fluctuations that accumulate along the propagation path, 
z, k is the wavenumber, and C,"(z) characterizes the magnitude of the optical 


turbulence (Fried, 1966). We now see the utility of the structure function approach in 


providing a robust parameter that can be used to characterize the atmospheric turbulence. 


2. Isoplanatic Angle 


Adaptive optical systems function by correcting the distorted wavefront based 
upon a measurement of the phase of the electromagnetic field. Since atmospheric 
turbulence varies with position, the correction for a specific source along an optical path, 


such as a star, is path dependent. The isoplanatic angle, 


L -3/5 
6, = a9 fc." @e"*ae (26) 
0 . 


is a measure of the angular separation between two sources for which the accumulated 


phase distortions are correlated (Fried, 1982). It is the angle where the Strehl ratio of an 


' of perfect performance. The Strehl ratio is the 


ideal adaptive optics system drops to e7 
ratio of actual peak intensity to the diffraction limited peak intensity. Isoplanatic angle is 


calculated by a vertical integration through the atmosphere weighted by the height, and is 
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basically the spherical path coherence length divided by the path length, 0, ~r,/L. The 


weighting occurs because two points higher in the atmosphere are farther apart and 


accumulate larger phase differences between the two propagation paths. 


35 Greenwood Frequency 
Wind advects the turbulent air across the optical path, changing the turbulent 


fluctuations. The Greenwood frequency, 


i= onoast? oy." (27) 


is found by integration over the optical path, taking into account the displacement caused 
by the wind component, v, , perpendicular to the integration path (Greenwood, 1977). It 
is a measure of the sampling rate an adaptive optics system must have to compensate for 


the turbulent degradation. Typically, an adaptive optics system with a bandwidth of 4 f 5 


will have a 3 dB loss in Strehl ratio. 


LY 
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Til, IMPLEMENTATION 


A. BALLOON FLIGHT 


A meteorological balloon sounding was taken at 0700 UT on 5 March 1995 from 
Wichita, Kansas, to measure the vertical temperature profile with high resolution (D. L. 
Walters, NPS). High-resolution and medium-resolution temperature instruments were 
added to the standard VIZ balloon package, providing 1.54 cm and 5.82 cm sella 
resolution temperature profiles, respectively. The theory of structure functions allowed 
subsequent computation of the optical turbulence profiles. 

The weather consisted of light, dry, snow below 2 km, topped by clear air. The 
balloon ascended through the stratosphere at 4.7 m/s on average, to an altitude of over 19 
km. Pressure, temperature, and humidity were recorded every 1.25 seconds, at an 
average 5.8 m spacing. Height above ground was determined from the hypsometric 
equation (Holton, 1979, pp. 19). Wind speed and direction were recorded every 30 
seconds. 

The standard VIZ low-resolution temperature probe had a time constant of about 
30 seconds and was sampled every 1.25 seconds. A medium-resolution, 100 jim bead 
thermistor added to a spare channel in the VIZ package provided more detailed 
temperature features of the atmosphere. The time constant increased with decreasing 


atmospheric pressure from 30 msec at ground level to 100 msec at 15 km. Hardware 
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limitations in the VIZ package limited the sampling interval to 1.25 sec. The temperature 
resolution of the micro-bead thermistor was between 0.005 and 0.01 K. (Fuehrer, 1994) 
Figure 5 shows the low-resolution and medium-resolution temperature profiles. 
The mismatch was caused by a non-linearity in the medium-resolution probe calibration, 
resulting in a negligible ~3% error. There was an inversion layer at 1.2 km and the 
tropopause was at 9.5 km. Figure 6 shows the medium-resolution potential temperature, 
which is the Kelvin temperature corrected for a dry adiabatic lapse rate, with 1000 mbar 
reference. The steep positive slope above the tropopause makes the stratosphere 
buoyantly stable. Large steps are visible throughout the profile. The wind speed and 
azimuth are plotted in Figure 7. High winds extended from 6 to 15 km, with strongest 


wind shear at the top boundary. 


B. HIGH RESOLUTION INSTRUMENT 


1. Description 

The high-resolution temperature probe was a 12.44:m_ chromel-constantan 
thermocouple, sampled approximately 315 times per second. Frequency modulation of 
the VIZ radio signal allowed the high data rate. This resulted in an average height 
resolution of 1.54 cm. The time constant of the thermocouple increased with decreasing 
pressure from 5 msec at ground to 15 msec at 20 km (Roper, 1992, pp. 62). To achieve 
milli-Kelvin resolution with a 12-bit A/D converter required a means of keeping the 


temperature in the middle of the dynamic range. This was accomplished by measuring 
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Figure 5. Low-resolution and medium-resolution temperature plotted against altitude. 
The mismatch between the profiles was caused by altitude non-linearity in the medium- 


resolution instrument. An inversion layer existed at 1.2 km and the tropopause was at 9.5 
km. | 
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Figure 6. Medium-resolution potential temperature plotted against altitude. The steep 
slope above 9.5 km makes the stratosphere buoyantly stable. The largest steps are visible 
throughout the profile. The troposphere was weakly stable, except for the strongly stable 
inversion layer from 1.2 to 1.9 km. 
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Figure 7. Wind speed and direction plotted with altitude. High winds extended from 6 to 
15 km, with strongest wind shear at the top boundary. 








the difference between the high-resolution thermocouple and a much slower responding 
thermocouple, with 10-30 second time constant from ground to 15 km. (Walters, 1995) 
Figure 8 shows the high-resolution temperature profile. Data was lost from 1.2 
km to 1.9 km because the A/D converter reached full scale in the sharply warming 
inversion layer. The stratosphere had larger temperature steps than the troposphere. 
Several regions ee large vertical temperature gradient were compared to confirm the 
conversion from A/D counts to Kelvin. Figure 9 is a plot of high and medium-resolution 
temperature for one of these regions. Both instruments agreed on the magnitude of the 


temperature change, confirming absolute calibration of the high-resolution instrument. 


2. Interference Spikes 


The 1.25 second interval read-out of low-resolution temperature by the VIZ 
package induced a signal into the co-located high-resolution cable. The resulting 
interference spikes are visible in Figure 9, spaced every 6 m, and were characterized by a 
sudden drop and rise in temperature lasting ~0.1 seconds. The magnitude of the spikes 
was about 0.02 K with a vertical extent of ~0.5 m. 

The effect of the spikes on the measured value of temperature turbulence was 
quite significant in regions of low turbulence. In a study of six non-turbulent regions 


throughout the stratosphere, ranging from 65 m to 185 m in vertical extent, a 3 to 10 
times increase in the local average C; due to the spikes was observed. The maximum 


effect occurred for sampling interval r of 30 cm to 40 cm, which was tuned to the width 
of the spikes. Significant contamination for sampling intervals of less than 2 m required 


their removal. 
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Figure 8. High-resolution temperature plotted against altitude. Data loss occurred from 
1.2 km to 1.9 km when the A/D converter reached maximum counts in the inversion 
layer. The stratosphere had larger temperature steps than the troposphere. 
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Figure 9. High-resolution temperature and medium-resolution temperature plotted 
against altitude. This temperature step was used to verify the conversion of high- 
resolution temperature from counts to Kelvin. The periodic interference spikes in high- 
resolution temperature are visible at approximately 6 m intervals. 
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The spike contribution to the temperature turbulence was removed by manually 
viewing the data and flagging the spikes. Numerous regions with turbulent fluctuations 
20.02 K obscured ~27% of the spikes and no attempt was made to remove them. Their 
contribution was considered negligible if not visually detectable. 

Since approximately 7% of the raw high-resolution data were discarded by spike 
removal, two processing procedures were used to compensate for the missing data. For a 


given sampling interval, +, four profiles of C7 were generated, each one offset by a 


quarter of the sampling sisal Subsequent averaging of the four C? profiles greatly 
reduced the occurrence of missing data. Figure 10 shows the histogram of the average 
value of C; in the stratosphere determined from 0.5 m sampling, where noise 
suppression was applied to the raw profile per Section 3. The distribution is visually log- 
normal. Based upon this observation, the remaining missing data were filled in with a 
local logarithmic average of 11 values. This number was adequate for the smallest r, yet 
not too large to alter the distribution significantly. Replacing the missing values 
maintained the approximately log-normal distribution. The procedure narrowed the 
distribution, reducing the linear average by 11%. 

The average value of Cc. in the stratosphere was computed for r less than 3 m to 
investigate the effect of the interference spikes on the stratospheric average. Figure 11 
shows that the spikes contributed to the turbulence mainly on length scales below 2 m. 


This confirmed again that the maximum effect was for 30 cm to 40 cm. 
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Figure 10. Histogram of the turbulence parameter in the stratosphere using a 0.5 m 
sampling interval. The effect of replacing data in gaps created by the spike removal is 
shown. Replacing the interference spike regions with a local logarithmic average 
narrowed the distribution and reduced the linear average by 11%. 
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Figure 11. Average turbulence parameter in the stratosphere plotted against sampling 
interval. Removing the interference spikes and replacing them with a local logarithmic 
average reduced the contamination, which was primarily below 2 m. 
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3. Noise Suppression 

The power spectral density of the high-resolution temperature in the stratosphere, 
with linear trend removed, is plotted in Figure 12. The slope of the power spectrum is 
between —2 and —3, and decidedly steeper than the —5/3 Kolmogorov result. The aideate 
spectral features from 1-10 Hz are harmonics of the 1.25 second period interference 
spikes. The sensor noise floor dominated above 30 Hz. The thermocouple time constant 
had a 3 dB cutoff at ~67 Hz, and the 315 Hz sample rate ensured oversampling of the 
probe. The power rose above 80 Hz from aliasing of the system white noise. 

To suppress noise, the raw data set was convolved with a triangular filter. This 
was equivalent to applying a sinc” filter in the frequency domain, with 3 dB point at 45 
Hz and first zero at 105 Hz. Noise was greatly reduced with the trade off that the 
minimum sampling interval, r, was limited to the 6 data point width of the triangle, or 
9.24 cm. 

The average value of the second-order structure function computed at the raw data 
resolution served as the noise estimate. Following the theory in Chapter II, Section C, the 


temperature can be written as a function of height, z , and time, r, 


T(z,t) = E(z,t)+ N(z), (28) 
where N(z) is the noise signal and &(z,t) the temperature turbulence. The presence of a 


slope discontinuity was ignored for simplicity. Assuming the system noise and turbulent 


fluctuation terms were independent leads to an equation similar to Equation (22), 
Dr (r)=([N, -2N, + NoP)+([E, -28, +4). (29) 
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Figure 12. Power spectral density of the stratosphere plotted against frequency. The 
power spectrum was steeper than the —5/3 Kolmogorov theory, shown with the dashed 
line. The distinct features from 1-10 Hz are harmonics of the 1.25 second period 
interference spikes. The noise floor dominated above 30 Hz and the power rose above 80 
Hz from aliasing of under-resolved system noise. 
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The noise is independent of the sampling interval, r , and so may be used for all r once 
determined. Because of the high-resolution instrument’s time constant, it was assumed 
that negligible real turbulent signal existed at the 315 Hz sampling rate. The value of the 


second-order structure function at 315 Hz was thus used as the noise estimate for all r. 
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IV. RESULTS 


A. TEMPERATURE STEP EFFECTS 


Temperature steps, or slope discontinuities, were abundant in the high-resolution 
temperature profile, recall the example in Figure 9. They also occur at the boundaries of 
“sheets” that have been observed by others, (Bertin, 1997), (Dalaudier, 1994), and 


mechanisms of their formation have been proposed (Coulman, 1995). Equation (22) in 


Chapter II describes their contribution to the second-order structure function, D); (r). 


The temperature structure constant, C,”, increased with sampling interval, r , see Figure 
11. This result confirmed that the Kolmogorov relationship for the structure function in 
Equation (7) failed for this data set. The positive slope of Ge versus r also indicates a 
power law with exponent greater than 2/3. 

Figure 13 shows the second-order structure function, Dr (r), calculated by 
averaging the entire stratosphere, where noise suppression and spike removal were 
performed per Chapter III. Both the high-resolution and medium-resolution data are 
plotted. Beyond 100 m, the A/D converter high-pass difference filter reduced the high- 
resolution signal. The Kolmogorov 2/3 slope and the step-removed values are plotted for 
reference and are further discussed in Section 2. 

The power law exponent, m(r), of D (r) is plotted in Figure 14. The high- 
| resolution and medium-resolution instruments agreed well below 100 m, the seh nite 


useful value of the high-resolution instrument. The rise in slope towards 2 for small 7 is 


35 








High—Resolution 

Hi-Res, Steps Removed 
10° Medium—Resolution 

Med—Res, Steps Removed 

Kolmogorov Theory 


Average D, (K’) 
o 
N 





10. 109 10! 10° 10> 
r (m) 


Figure 13. Average 2" order structure function of the high-resolution temperature profile 
in the stratosphere (10.4-19.4 km) plotted against sampling interval. Results of step 


removal follow the Kolmogorov theoretical 2/3 slope until ~10 m. High and medium- 
resolution data match well. 
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Figure 14. Slope of the 2" order structure function plotted against sampling interval. 
High and medium-resolution slopes match well below 100 m. 
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evidence of the inner scale length, /,, coupled with the onset of the viscous convective 


region. Note that the 3 dB point of the thermocouple was at ~4.7 cm. The effects of 
dissipation were quite significant below ~25 cm, a value for which the slope has risen 
halfway from the flat region of the curve to the theoretical value of 2. This marks the 
inner scale of the inertial subrange and corresponds to ~20 Hz on the power spectrum in 
Figure 12. Averaging the stratosphere smears the inner scale effect because the inner 
scale grows with decreasing pressure at higher altitudes. 

The slope falls rapidly with r until reaching a plateau from 0.5 mto2m. A 
steady rise in slope above 1 m indicates the contribution of step discontinuities to the 
structure function. Apparently, the steps prevented the slope from ever reaching 2/3. The 
Q ratio of Equation (25) predicted an increase in the effective exponent, mr). It was 
expected that the slope would rise to about 2 and level off when steps produced the 
majority of the structure. The fact that it did not do so implies the step effect was 
reduced above 40 m. This behavior may characterize the maximum size of small steps. 
The rise in slope beyond 200 m possibly indicates a regime of large steps. The analysis is 
limited for large r because the number of values used in the average declines inversely 


with r. 


B. REMOVAL OF STEPS 


From the behavior of D");(r), it was initially concluded there was no 


Kolmogorov inertial subrange along a vertical profile of the stratosphere. Based upon the 
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theory of steps presented in Chapter II, an attempt was made to eliminate uke 
contribution to the structure function. The 2" and higher order structure functions 
mathematically fail in the presence of rapid slope pines or temperature steps, leading 
to a criterion for their identification. Visual judgement by the investigator determined the 
presence of reasonably obvious slope changes that could not be random turbulence. The 
rapid slope change across the step temperature increase in Figure 9 is an example at 5 m 
resolution. At 100 m resolution this particular step was not resolved and did not cause a 
significant shift in the slope. The change in slope was very gradual at the minimum 
useable data resolution of 9.24 cm, therefore this step was not removed. The main 
vertical extent of the turbulent sheets in the stratosphere ranges from meters to tens of 
meters (Dalaudier, 1994). Only the top and bottom edges of these sheets were deemed 
steps and removed by this procedure, leaving the middle of the sheet intact when the 

sampling interval permitted. 
_ The steps were identified for several selected sampling intervals for both 


instruments, listed in Table 1. The investigator manually viewed every data point 


High-Resolution Medium-Resolution 
wm 


| 0.944e-4 


0.903e-4 


l.1lle-4 
1.39e-4 
1.67e-4 





Table 1. Results of step removal. 
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at each of the selected r values and flagged the affected regions. All flagged data were 


replaced with the local logarithmic average in the manner of processing the interference 
spikes in Chapter III. Uncertainty in C,’ is reduced by the 3/5 power when calculating 
the optical parameters 7,, 6,, and f 2: Complexity of developing an automated 
algorithm along with limited time necessitated a manual procedure. The Q value could 
prove useful in developing an automated routine by defining a maximum permitted value, 


Ona» he maximum allowable slope difference would then be, 


B. =|2.41C29. | 23, (31) 
max T max 
Figure 15 is a comparison of the ce profiles after step removal for the 5.08 m 


high-resolution and 5.82 m medium-resolution cases, with 250 m smoothing for 
visualization. The profiles agree well except above 18 km, where both probes produced a 
similar turbulence spike at 18.5 km, although offset by roughly a factor of 3. Excessive 
manual step removal in the medium-resolution profile was most likely the cause of this 


discrepancy. In addition, there is a good correlation between high turbulence and the 
strong wind shear regions in Figure 7. The pre and post step-removal dependence of C - 


on r is shown in Figure 16. The two leftmost step values were averaged to obtain 


(C,’)=0.924e-4 (K*m’*'’). These two points were chosen because they appeared to 


represent an asymptotic limit of the step removal process. It is interesting that this value 
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Figure 15. Turbulence parameter in the stratosphere plotted against altitude. The steps 
were removed for sampling interval r = 5.08 m for high-resolution and r = 5.82 m for 
medium-resolution. High and medium-resolution profiles match well below 18 km. A 
250 m smoothing filter was applied for visualization. 
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Figure 16. Average turbulence parameter in the stratosphere (10.4-19.4 km) plotted 
against sampling interval. Manual removal of steps was performed for 11 sampling 
intervals. The curve formed by these points has an asymptotic value of 9.24e-5 for the 
turbulence parameter. The approximately constant region corresponds to a Kolmogorov 
inertial subrange with outer length scale L,~10 m. 
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intersects the unprocessed curve at ~25 cm. The step-removed values depart from a near 


constant value above 10 m, making this an estimate of the outer scale length, L,, of the 


Kolmogorov inertial subrange. The step removed D”), (r) values plotted in Figure 13 
also depart noticeable from the 2/3 theoretical slope above 10 m. 

The four bends in Figures 13, 14, and 16 summarize the physics of the 
atmosphere for this balloon launch. The bend at 25 cm marks the inner scale boundary 
between the dissipation regime and the inertial subrange. A bend at 3 m marks the scale 
where temperature steps become the dominant source of structure and steadily grow in 
importance, according to the Q ratio in Equation (24). Steps came into play at small r, 
before a Kolmogorov inertial subrange could be established. The decrease in slope at 40 
m characterizes a scale length above which step effects declined. A final increase in 
slope at 200 m indicates an increase in the steps, although the statistics are not as good 
for this range of sampling intervals. This could indicate that the distribution of step sizes 
has multiple modes, with a gap between 40 m and 200 m. 


Removal of steps revealed an underlying outer scale length, Z,~10 m. Below L,, 
the procedure uncovered a nearly Kolmogorov inertial subrange. Earlier research noted 
that Cc. from a high-resolution vertical profile, r=1 m, was a factor of ~1.6 greater than 


that measured by a 1 m horizontal pair of temperature probes (Walters, 1995). Figure 16 
predicts a factor of ~2. Since horizontally spaced probes are also susceptible to steps, 
from pendulum motion of the package under the balloon in the presence of vertical slopes 


in the stratified layer, this investigation indicates 1 m horizontally spaced probes over- 
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predict turbulence by perhaps 10-25%. The inability of horizontal probes to measure a 
horizontal temperature profile makes the step removal process impossible with current 
systems. The application of a correction factor is unwarranted without further research to 
quantify the factor more accurately. 

The results are based upon the stratosphere only, but the theory should also be 
applicable to the troposphere. Current meteorological weather balloons, with the addition 
of an inexpensive bead thermistor, would yield a temperature profile of ~5 m vertical 
resolution, though 1-2 m would be preferred. Under lower turbulence conditions than 
this investigation, the stratified temperature profile would be less disturbed and the step 
contribution is expected to be lower; the procedure should continue to work. Higher 
turbulence is expected to raise the outer scale length, while simultaneously increasing the 
effect of steps. Further research should reveal which of these two competing effects 
dominates and show if a 5 m vertical resolution temperature profile is sufficient for a 


wide range of atmospheric conditions. 
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Vv. CONCLUSIONS 


Spatial variations in the index of refraction of the atmosphere distort propagating 
electromagnetic waves. Many laser and telescope systems of interest to the military 
operate within the visible through near infrared portion of the spectrum and are 
susceptible to these effects. The effects of atmospheric turbulence are the primary 
limitation to performance of high angular resolution systems. The Air Force’s Airborne 
Laser program and other military adaptive optics systems would benefit from the ability 
to forecast severe atmospheric turbulence situations accurately. Good validation data is 
required to develop such a model. 


Structure functions provide a robust means of extracting the index of refraction 


2 


n ? 


structure constant, C,°, a useful turbulence parameter in the atmosphere. Vertical 


temperature profiles through the atmosphere reveal the presence of steps in the potential 


temperature. The rapid change in slope of the temperature across these steps contributes 
a false term to C,’. This research developed a quantitative measure of the false 


turbulence contribution caused by the temperature slope discontinuities in the 
atmosphere. The homogeneity and isotropy requirements of structure functions led to a 
technique of removing the false contribution. 

Both 1.54 cm and 5.82 m vertical resolution profiles of temperature were 
collected from a balloon flight in March 1995 from Wichita, Kansas. Data processing 
removed and compensated for interference spikes in the high-resolution temperature 


profile. Treating the stratosphere as an ensemble average, the application of a second- 
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order structure function to the temperature profiles confirmed the presence of a 
temperature step contribution. Removal of the temperature steps for several sampling 
intervals exposed a Kolmogorov inertial subrange extending from ~25 cm to ~10 m. 

The results were based upon the stratosphere only, but the theory should be 
applicable to the troposphere, where the temperature steps are less pronounced. Future 
research could extend the analysis to the troposphere. Current meteorological weather 
balloons, with slight modification, are capable of measuring temperature profiles with 5 
m vertical resolution. Based upon the results of others, wind shear appears to be the 
source of the steps (Bertin, 1997) (Dalaudier, 1994). The effect of these steps should be 
less pronounced under lower turbulence conditions, and the procedure should continue to 
work. Higher turbulence is expected to raise the outer scale length while simultaneously 
increasing the magnitude or number of steps. Further research should decide if 5 m 
véttival resolution is sufficient for a wide range of turbulence conditions. 

This investigation produced a reference turbulence profile in the stratosphere for 
validation of turbulence forecasting models. The potential of this technique to compute 


the isoplanatic angle, 8, , coherence length, r,, and Greenwood frequency, f,; reliably 


with inexpensive balloons should benefit ABL and other national programs. 
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